Mechanism behind self-sustained oscillations in direct current glow discharges and 

dusty plasmas 



o 

CN 

Or 
< 



43 



c*3 
in 

O 

• i-H , 

>> 
43 ■ 



CN 
> 

o 
l> 

CN 
O 
cn 



Sung Nae ChcQ 

Devices R&D Center, Samsung Advanced Institute of Technology, 

Samsung Electronics Co., Ltd, Mt. 14-1 Nongseo-dong, 
Giheung-gu, Yongin-si, Gyeonggi-do 446-712, Republic of Korea. 
(Dated: 2 April 2013) 

An alternative explanation to the mechanism behind self-sustained oscillations of ions in direct current (DC) 
glow discharges is provided. Such description is distinguished from the one provided by the fluid models, where 
oscillations are attributed to the positive feedback mechanism associated with photoionization of particles and 
photoemission of electrons from the cathode. Here, oscillations arise as consequence of interaction between 
an ion and the surface charges induced by it at the bounding electrodes. Such mechanism provides an elegant 
explanation to why self-sustained oscillations occur only in the negative resistance region of the voltage-current 
characteristic curve in the DC glow discharges. Furthermore, this alternative description provides an elegant 
explanation to the formation of plasma fireballs in the laboratory plasma. It has been found that oscillation 
frequencies increase with ion's surface charge density, but at the rate which is significantly slower than it does 
with the electric field. The presented mechanism also describes self-sustained oscillations of ions in dusty 
plasmas, which demonstrates that self-sustained oscillations in dusty plasmas and DC glow discharges involve 
common physical processes. 



I. INTRODUCTION 

When an ion is confined between the two plane-parallel 
electrodes and is subject to static electric field (see Fig. 
[T]), it begins to self-oscillate without damping. Such phe- 
nomenon is referred to as self-sustained oscillations in di- 
rect current (DC) glow discharges; and, it has been known 
for almost a century^"— In literature, the phenomenon of 
DC glow discharge is also referred to as the DC glow 
corona. The mechanism behind such oscillations is still not 
fully understood] 8 ' 9 ! 1 l ~ J2. Over the years, various theoreti- 
cal models have been proposed in an attempt to explain the 
phenomenon ! 5 ' 1 M 8 " — Among the successful ones are those 
based on the fluid and equivalent circuit models. 

The equivalent circuit models try to predict oscillations by 
representing the system with an equivalent RLC circuit, where 
R is a resistor, L is an inductor, and C is a capacitofr!^"— Al- 
though this approach is quite useful in describing oscillation 
frequencies as function of DC bias voltages, it says nothing 
about the mechanism behind self-sustained oscillations. 

The fluid models approach the problem from more fun- 
damental grounds of the electromagnetic theory. In this ap- 
proach, the Poisson equation is solved in combination with 
the electron and the ion flux continuity equations, which con- 
stitute the so called positive feedback mechanism. 22 " 26 Asser- 
tion of such feedback mechanism is crucial in the fluid models 
because, without it, no oscillatory solutions can be obtained. 
Typical sources of the positive feedback mechanism include 
the photoionization of particles and photoemission of elec- 
trons from the cathode. Physically, such feedback mechanism 
promotes oscillations by periodically reversing the particle's 
charge polarity. An ion oscillating near an electrode induces 
current in the same electrode, where the waveform of such 
current is correlated to its velocity profile. 23 Experimentally, 
it is this induced electrode current (or voltage) which gets 
measuredi&ii The hallmark of the fluid models is their poten- 
tial to reproduce experimental measurements to a reasonably 



good accuracy. In principle, with an appropriate specification 
of the positive feedback mechanism, the level of accuracy pre- 
sented by the fluid models can always be improved. The feed- 
back mechanism varies among different authors; and, this has 
been a subject area of ongoing debate among different fluid 
model theorists ! n ' 22 ~ — 




Figure 1 : (Color online) Illustration of core-shell structured ion con- 
fined between DC voltage biased plane-parallel conductors with a 
vacuum gap of h. E p = — e z (Vr — Vl) /h is the parallel plate electric 
field, where (e X) e y ,e z ) are versors along the Cartesian (x,y,z) axes, 
respectively. 

Does the aforementioned positive feedback mechanism, 
which is asserted in the fluid models, represent the funda- 
mental mechanism behind the phenomenon of self-sustained 
oscillations in the DC glow discharges? This is a subtle ques- 
tion because I have shown recently that an ion confined be- 
tween the DC voltage biased plane-parallel conductors goes 
through an undamped self-oscillatory motion. 27 Such self- 
oscillatory motion requires that an ion is electrically polar- 
izable, but it does not necessarily involve or require the dis- 
cussed positive feedback mechanism which is asserted in the 
fluid models. The fact that an ion must be electrically po- 
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larized excludes electrons from consideration in the discus- 
sion of self-oscillations presented in this work. However, the 
proton, which is known to be electrically polarizable, 28 is ex- 
pected to self-oscillate in the DC glow discharges according 
to this alternative theory. 

Remarkably, the predictions of this alternative theory^ 9 , 
qualitatively agree with the predictions made by the fluid 
models. 23,25 Both theories predict a saw-tooth shaped wave- 
form for the induced currents in the electrodes. Furthermore, 
sharp pulses of radiation output are predicted by both theories 
to accompany the abrupt rises in the induced electrode cur- 
rents. Such remarkable similarities in the predictions by the 
two very dissimilar theories suggest that the discussed positive 
feedback mechanism, which is asserted in the fluid models, 
may not necessarily represent the fundamental mechanism be- 
hind the phenomenon of self-sustained oscillations in the DC 
glow discharges. 

In this work, I shall present some new aspects of this al- 
ternative theory 29 for further exploitation of the phenomenon 
of self-sustained ion oscillations in the DC glow discharges. 
These theoretical predictions are explicitly compared with 
various experimental results for the validation of the model. 
By direct application of the model to describe the charged- 
particle oscillations in dusty (or complex) plasmas experi- 
ments, I shall reveal that the phenomenon of self-sustained 
oscillations in both dusty plasmas and DC glow discharges 
involves a common physical mechanism in which self- 
oscillations are attributed to the interaction between an ion- 
ized particle and the surface charges induced by it at the sur- 
faces of the particle confining electrodes. Considering that 
particle oscillations in both dusty plasmas and DC glow dis- 
charges experiments involve ionized particles that are only 
differ in their physical sizes and masses, such outcome is not 
too surprising. After all, apart from this, everything else is 
nearly identical from the physics point of view in both exper- 
iments. 



tive theory using the configuration illustrated in Fig. Q] The 
electric potential between the plates is obtained by solving 
the Laplace equation with appropriate boundary conditions. 27 
With this electric potential, induced surface charges at the sur- 
faces of conductor plates are obtained by application of the 
Gauss's law. Such charge distributions act on an ion and gen- 
erate resultant force given by F7- = F] + F2 — e z mg, where m 
is the mass of an ion, g = 9.8 m/s 2 is the gravitational con- 
stant, and F] (F2) is the force between the ion and the surface 
charges induced by it at the surface of anode (cathode). The 
explicit forms of Fi and F2 are given by 2 - 7 - 
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where Co is the vacuum permittivity, E p = ||E p || is the 
parallel-plate electric field strength, the parameter s is the dis- 
tance from the ion's physical center to the anode; and, the 
terms y and v are defined as 
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II. BRIEF OUTLINE OF THE THEORY 

The self-sustained oscillations in DC glow discharges are 
now discussed briefly in the framework of presented alterna- 



where <Ti (02) is the surface charge density at the ion's core 
(shell) of radius r = a (r — b), and the dielectric constants K2 
and K3 are depicted in Fig. Q] The resultant force on an ion is, 
hence, given by 
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where s = id + b in Fi and F2 (see Fig. [TJ. The potential energy associated with this force is given by 29 
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and, the equation of motion associated with this force, i.e., Eq. (0), can be expressed as 
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where /and v are defined in Eqs. ((TJ and (fJJ, respectively. 27 ' 29 

III. COMPARISON TO THE EXPERIMENTS 

Gases under pressure in general, including noble gases, 
are composed of atomic clusters, a state of matter between 
molecules and solids, due to Van der Waals interaction ! 30 " 32 
In gaseous argon, each spherical atomic clusters of radius 
r w lnm contains roughly 135 argon atoms. 31 Modeling par- 
ticles in the gas as atomic clusters is important because gases 
are delivered to the laboratories in pressurized containers, in 
which environment individual particles in gas exist in the form 
of atomic clusters. That said, I shall explicitly work with an 
argon atomic cluster of radius r « 1 nm, which for brevity 
is simply referred to as "ion" hereafter. Such ion is not ex- 
pected to be core-shell structured like the one depicted in 
Fig. 03 The shell portion of an ion can be eliminated by 
choosing b = a, K2 — °°, and 02 = OC/m 2 . For a positive 
ion, its surface charge density is given by <J\ = Nq e / (4-na 2 ^ , 
where N is the number of electrons removed from the parti- 
cle and q e = 1.602 x 10~ 19 C is the fundamental charge mag- 
nitude. Without loss of generality, and for the purpose of 
clear illustration in this work, I shall choose N = 250 and 
a = 1 nm. This corresponds to the ion's surface charge den- 
sity of <7i ss 3.19C/m 2 . Assuming the mass of an argon atom 
is 6.63 x 10~ 26 kg, a spherical atomic cluster composed of 135 
argon atoms has a total mass of m w 8.95 x 10~ 24 kg, where 
the masses of missing N electrons have been neglected. Purely 
for convenience, I shall assume that the cathode is grounded 
and the space between the two conductor plates in Fig. [T] is 
a vacuum with a gap of h = lOOnm. The obtained parameter 
values are summarized here for reference: 

K 2 =°°, JC 3 = 1, /!=100nm, V L = 0V, 

cti «3.19C/m 2 , a 2 =0C/m 2 , (6) 
m k, 8.95 x 10~ 24 kg, b = a = lnm. 

Illustrated in Fig. [2] is a plot of the potential energy well, 
where Eq. has been plotted for Vt = 100 V using the 
parameters defined in Eq. ©. The width Id of the poten- 
tial energy well decreases with an increase in the electric field 
strength E p . Such characteristic provides an elegant explana- 
tion to the experimental observations in which the oscillation 
frequencies increase with the applied electric field strength. 
Fox seems to be the first to discuss on such characteristic us- 
ing a parabolic potential model.— 

In the case of positive ions, the potential energy well gets 
formed in vicinity of the anode at the presence of applied 
static electric field. On the other hand, for the negative ions, 
the same electric field results in formation of the potential 
energy well in vicinity of the cathode. Such characteristic 
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Figure 2: Potential energy is plotted for Vt = 100 V using parameter 
values defined in Eq. @. Here, Vj = 100V corresponds to E p = 
1 GV/m. 



provides an elegant explanation to the phenomenon of fire- 
ball formation in the laboratory plasma. When a positively 
biased electrode is immersed in plasma, a region of intense 
glow appears in vicinity of the electrode. Such glow forma- 
tion is referred to as a plasma fireball. 33 In the framework of 
this alternative theory 27 ' 29 elaborated in this work, the forma- 
tion of such plasma fireballs can be elegantly explained from 
the potential energy well which gets formed near the anode 
for positive ions. Since the dynamics of a plasma involves 
collective motions of its constituent atoms, a weakly ionized 
plasma fireball, as a whole, self-oscillates near the electrode 
as if it were a weakly ionized single superparticle. Through- 
out this work, I shall use the term "superparticle" to refer to 
such entity as plasma fireball whose dynamics can be repre- 
sented by an equivalent single-particle picture. In the equiva- 
lent single-particle picture, the self-oscillation dynamics of a 
plasma fireball is described by the equation of motion defined 
in Eq. © with appropriate effective mass and surface charge 
densities prescribed. With that in mind, the plasma fireball ef- 
fectively has a very large mass and carries a net charge which 
is very weak. Consequently, induced current oscillations of 
relatively low frequencies and small amplitudes are expected 
to be generated at the anode near a self-oscillating plasma fire- 
ball. Such prediction is consistent with the experimental ob- 
servation by Stenzel et a/. 33 Notice that this alternative theory 
can be directly applied to describe the self-oscillations of ion- 
ized dust particles in dusty plasmas. It is remarkable that the 
profile of the potential energy well illustrated in Fig. |2]quali- 
tatively agrees with the measurements by Tomme et al^l and 
Arnas et alM in dusty plasmas. Further discussion on this is 
provided in the sections A thru C of the appendix. 

The particle's equation of motion described in Eq. ((5} can 
be solved via Runge-Kutta method using the parameter values 
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defined in Eq. (O and the initial conditions given by 

z rf (0)=0.5nm and ^=0^. (7) 

The results for Vj — 50 V, 75 V, and 100 V are shown in 
Fig. Oa), where it shows the oscillation frequencies increas- 
ing with the electric field strength. The oscillation frequen- 
cies also increase with the ion's surface charge density, but 
at the rate which is significantly slower than it does with the 
electric field strength (see Fig. 0b)). For instance, when 
the electric field strength is doubled from E p = 0.5GV/m 
to E p = 1 GV / m, oscillation frequency nearly doubles from 
V sc ~ 104 GHz to v osc « 191 GHz. However, doubling the 
ion's surface charge density, i.e., (7i = Nq e / (4na 2 ) , from 
N = 250 to N = 500 only slightly increases the oscillation fre- 
quency from Vosc ~ 191 GHz to v OS c ~ 207 GHz. Such charac- 
teristics is consistent with the observations by Fox 3 and Bosan 
et aJ., 36 where they have reported that oscillation frequencies 
do not seem to be very dependent on the type of ions used in 
the discharge. It is well known that dissimilar atomic gases 
have different ionization tendencies. Based on this, it can be 
inferred that Ci , in general, for different atoms under identical 
conditions are different. 

Although a relatively high static electric field is being ap- 
plied to the ion in Fig. EJa), its kinetic energies are small 
due to the fact that ion is going through frequent turnings in- 
side a tight potential energy well. The ion's kinetic energy 
decreases with increased electric field strength for that matter 
and this is illustrated in Fig. EJc). At the maximum speed of 
~ 2.912km/s, the ion gains kinetic energy of s» 237 eV, 
which is small compared to the depth of pt, rs 4.1 keV for the 
potential energy well but large compared to the average room 
temperature thermal energy of Sj = l.SkgT « 39meV (here, 
T = 300K and k B = 8.62 x 10~ 5 eV/K is the Boltzmann con- 
stant). Self-sustained oscillations in the DC glow discharges 
can thus be expected to have high thermal stability. 

In the experiments, self-sustained oscillations of ions in the 
DC glow discharges occur only in the negative resistance re- 
gion of the voltage-current characteristic curve . 12,14 ' 18 ' 24 Such 
negative resistance region in the voltage-current characteris- 
tic curve is characterized by Vthi < Vr < V"th2, where Vthl 
and V t h2 are some threshold voltages. Equation (]5) provides 
an elegant explanation to such properties. For instance, Eq. 
(O yields solutions that are unphysical for E p < Ep^\ and 
E p > E p fb2, where En,thi ar, d £p,th2 are the threshold elec- 
tric field strengths. The physical solutions are only obtained 
for Zip thi < E p < Epfia- It is difficult to pinpoint E pX hi and 
Ep t h2 without the analytical solution of Eq. ((5) at hand. Nev- 
ertheless, these can be roughly estimated on the grounds of 
physicality. To illustrate this, the oscillation frequency is first 
plotted as a function of the electric field strength in Fig. |3d) 
by numerically solving Eq. (0 using the parameter values and 
the initial conditions from Eqs. (|6]l and Q, respectively. The 
solutions obtained for E p < O.lGV/m and E p > 4.9GV/m 
are unphysical. For instance, in the case of E p > 4.9GV/m, 
solutions show that the ion penetrates into the anode's surface 
whereas fovE p < 0.1 GV/m, solutions yield peak to peak am- 
plitudes that are larger than h/2. This latter case, although 



mathematically allowed, is unphysical because, for a posi- 
tive ion, oscillations can only exist for Zd < h/2 — b from 
the argument based on the grounds of physicality. — Hence, 
the physical oscillatory solutions exist only for 0.1 GV/m < 
e p S 4.9 GV/m, where the upper and the lower bounds of Ep 
are only rough estimates based on the grounds of physicality 
arguments. Further discussion on the properties of the phys- 
ical and the unphysical oscillatory solutions are provided in 
the section D of the appendix. 

Recently, Lotze et aJ. 37 reported on an indirect account of 
oscillations involving a single H2 molecule in the DC volt- 
age biased conductors near an absolute zero temperature of 
T = 5 K. According to their findings, the H2 molecule in the 
junction, which is the space between the atomically clean 
Cu(lll) surface and the STM (scanning tunneling micro- 
scope) tip mounted on a cantilever, self-oscillates between 
the two unknown positional states when a threshold DC bias 
voltage is applied to the electrodes, i.e., the copper Cu(lll) 
surface and the STM tip. The footprint of self-sustained os- 
cillations is the presence of the negative differential conduc- 
tance in their measurement. Gupta et al— showed that such 
negative differential conductance is due to the H2 molecules 
in the junction. In the DC glow discharges, self-sustained 
oscillations arise as a consequence of the negative differen- 
tial resistance . 1Q ' 12,14 Since the negative differential conduc- 
tance and the negative differential resistance are reciprocally 
related, they represent an equivalent description of the same 
physical processes. It is well known that the H2 molecules 
can be ionized by a process of an electron impact. 39 When 
a threshold DC bias voltage is applied to the electrodes, the 
energetic electrons begin to tunnel through the junction. It 
is possible that the H2 molecules in the junction are ionized 
in the process. If that is indeed the case, what Lotze et al. 
reported may be an indirect account of self-sustained oscilla- 
tions involving a single or a few ions in the DC glow corona. 29 
Their observation can serve as a catalyst for the future exper- 
iments in the few-particle DC plasma experiments. 

Arnas et al. 40 reported that a dust particle made of hol- 
low glass microsphere of radius r w 32 jj.m with mass density 
of p m ss 110kg/m 3 and carrying a total electrical charge of 
Q « -4.3 x \0 5 q e , where q e = 1.602 x 10~ 19 C, self-oscillates 
at a frequency of v osc ~ 17 Hz inside the plasma sheath. The 
plasma sheath environment, in many respects, is similar to 
the empty space region between the DC voltage biased plane- 
parallel conductors. Ionized particle oscillations in the plasma 
sheath can therefore be modeled from the simple apparatus il- 
lustrated in Fig. Q] That said, the parameters and the initial 
conditions in the experiment by Arnas et al. can be summa- 
rized as follow: 

K"2 = 6.5, K3 = 1 . /i = 5mm, E p = 5.15kV/m, 
cti =0C/m 2 , a 2 = -5.35 /iC/m 2 , 
m ~ 15.1 ng. a = 0m, b = 32jim, 
Zd(0)= 4.4mm, id (0) = Omm/s, i d = dz d /dt, 

where the details of how these were obtained are explained 
in the section B of the appendix. Using these, the equation 
of motion described in Eq. © can be solved numerically via 
Runge-Kutta method. The result is shown in Fig. [4] One 
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Figure 3: (Color online) Equation (f5} is plotted, where in (a) electric field is varied and in (b) ion's surface charge density d\ = Nq e / (47Ta 2 ) is 
varied at E p = 1 GV/m. (c) Ion velocity corresponding to the plot in (a), (d) Dependence of oscillation frequency v osc on electric field strength 
E p . Similar plot for the lower frequencies involving weakly charged ions subject to smaller E p is also provided in the section E of appendix. 
In (a,b,c,d), all of the unspecified parameter values are from Eq. l[6j and the initial conditions are from Eq. (0. 



can readily verify that the particle oscillates at a frequency 
of approximately v osc ~ 17 Hz, which is consistent with the 
measurement by Arnas et al£& 

Although Amas et al. explicitly states that their glass mi- 
crosphere is hollow structured, their paper does not provide 
any information regarding the radius of the void inside of 
their hollow glass microsphere. Due to the lack of this in- 
formation, the mass of an hollow glass microsphere in Eq. 
dHJ has been computed assuming a solid microsphere. Such 
assumption overestimates the mass of an hollow glass mi- 
crosphere. In fact, in their paper, Arnas et al. indicates the 
gravitational force of F g = (1.5 ±0.3) x 10~ 10 N for their hol- 
low glass microsphere. 40 Such force corresponds to a mass 
of m ~ 15.3 ng, which is exactly the mass computed as- 
suming a solid glass microsphere in Eq. (HJ. If instead a 
smaller mass of m m 6.19ng is assumed, an electric field 
strength of E p = 2.2kV/m is sufficient to generate an oscilla- 
tion frequency of v osc ~ 17 Hz. This value for the electric field 
strength, i.e., E p — 2.2kV/m, is consistent with the local elec- 
tric field strength in the plasma sheath obtained by Amas et al. 
in their dusty plasma experiment. The particle's oscillatory 
motion corresponding to this latter case, where m « 15.1ng 
and E p = 5.15kV/m, has been plotted in Fig. |4]for compari- 
son. 



m = 6.19 ng, E = 2.20 kV/m 
m = 15.1 ng, E, =5.15kV/m 
4.9 - ' 




0.05 0.1 0.15 0.2 0.25 0.3 
f[s] 

Figure 4: (Color online) Equation ((5} is plotted using parameter val- 
ues and initial conditions defined in Eqs. l[8}. Oscillation frequency 
of v osc ~ 17 Hz is predicted by the theory, which is consistent with 
the measurement by Arnas et al.^° In the plot, the anode is located at 
Zd = 0mm and the cathode is located at id — 5 mm. 



This remarkable result demonstrates that self-sustained os- 
cillations in both dusty plasmas and DC glow discharges share 
a common physical mechanism, in which the oscillations can 



6 



be attributed to the interaction between the ionized particle 
and the surface charges induced by it at the bounding elec- 
trodes. Such outcome makes sense because the oscillating 
ionized particles in both dusty plasmas and DC glow dis- 
charges experiments are only differ in their physical sizes and 
masses. Apart from this, the two systems are nearly identical 
in nature from the physics point of view. 

Quite often in dusty plasmas, theoretical models based on 
simple harmonic excitations of small amplitude oscillations 
are employed for the description of the particle dynamics. One 
such model is given by^- 

d 2 z dz ^ 2 /o cos co/ 

-TTI + Vdn-T- +£&=- , (9) 

dt L dt nii 

where z is the particle displacement, is the particle's mass, 
(0 is the angular frequency of oscillation, t is the time, and 
fo is the amplitude of the external force. Particle's equation 
of motion based on such model, Eq. (0, is limited to the 
descriptions of oscillations involving a small amplitude dis- 
placements about the stationary equilibrium position. Besides 
this limitation, the equation of motion described by Eq. (0, as 
it stands, cannot be directly applied to describe the particle's 
oscillatory motion due to the fact that the quantities Vd n and 
Q. v must be first determined experimentally. 42 Consequently, 
Eq. (0 is only useful when determining the total charge car- 
ried by the particle. Although Eq. (0 provides an indirect 
method of measuring the particle's total charge, as it stands, 
it says nothing about the fundamental mechanism behind the 
self-sustained oscillations of ions in dusty plasmas. 

The equation of motion described in Eq. (0 is distin- 
guished from the one illustrated in Eq. (0 in that its physi- 
cal description is not limited to just small amplitude oscilla- 
tions but covers particle oscillations of any amplitude ranges. 
Further, Eq. (0 does not depend on quantities like Vd n and 
Q. v , which terms must be determined experimentally. Instead, 
Eq. (0 completely describes the dynamics of charged-particle 
motion solely based on the information obtained from the 
particle's physical properties (i.e., mass, surface charge den- 
sity, and dielectric constant, etc.) and the local electric field 
strength. In the case where the mass of an ionized particle 
is represented by that of an ionized atom, Eq. (0 describes 
the dynamics of self-oscillating ionized atom in the DC glow 
discharges. In the previous work, 29 I have discussed that Eq. 
(0 yields results that qualitatively agree with certain aspects 
of self-oscillations in the DC glow corona experiments pre- 
dicted by theories based on the fluid models i 23 ! 25 If plasma 
in the DC glow corona experiment is effectively treated as a 
self-oscillating weakly ionized single superparticle, Eq. (0 
provides a satisfactory description of the electrode current os- 
cillations in the DC glow corona. Further discussion on this is 
provided in the section E of the appendix. 

Until now, no single theory was able to successfully ex- 
plain self-sustained oscillations in both dusty plasmas and DC 
glow discharges experiments. For years, experimental evi- 
dences hinted that these were related phenomena, but without 
any definitive conclusions. In this respect, Eq. (0 is the first 
theoretical model to provide such definitive conclusions. 



IV. DEVICE APPLICATION 

It is well known that oscillating ions generate electric dipole 
radiation. Such property can be utilized to develop a wideband 
electromagnetic radiation source in which the frequency of 
emitted radiation can be tuned by varying the DC bias voltage 
across the electrodes i 27,29 According to Zouche and Lefort, 
the plane-parallel plate electrodes made of nickel-silver com- 
posite material, which are separated by a gap of h = lOOnm, 
can support DC bias voltages up to ~ 400 V across the elec- 
trodes before an onset of electrical breakdown. 43 With im- 
provements in the electrical breakdown characteristics, such 
device can be engineered to cover both the microwave and 
the terahertz band of electromagnetic spectrum with high effi- 
ciency. 

The first experimental evidence that the phenomenon of 
self-sustained ion oscillations in the DC glow discharges has a 
potential applications in the development of a wideband elec- 
tromagnetic radiation source came from McClure in 1963, 
when he reported an oscillation frequency of v 0S c ~ 20 MHz in 
a low pressure glow discharge tube, which comprised of cylin- 
drical hollow cathode and a very thin coaxial wire anode. 44 
Unfortunately, McClure provided no theoretical model to ex- 
plain his finding. Such concentric cylinder configuration for 
the electrodes is a typical apparatus found in many DC glow 
discharges experiments. Nearly two decades later, in 1980, 
Alexeff and Dyer shrunk the aforementioned concentric cylin- 
der configuration (filled with air at gas pressure of 0.1 mTorr) 
to the size of a pen and successfully demonstrated the gen- 
eration of microwave radiations in the gigahertz frequency 
ranges. 45 Their pen sized coaxial configuration was later re- 
ported to generate radiation at terahertz frequency of V osc « 
1 THz with an output radiation power of P ra d ~ 1.5W.i^In 
order to explain their observations, Alexeff and Dyer propo sed 
a model which they referred to as the orbitron theory £2.~2lThe 
basic idea behind the orbitron theory is very simple. The elec- 
trons emitted from the inner surface of the outer concentric 
cylinder (which part represents the cathode in the coaxial con- 
figuration) orbit around a thin coaxial wire anode; and, such 
electron orbital motion generates the detected high frequency 
electromagnetic radiation. However, various experimental re- 
sults reported by other groups 48 " 50 contradicted the orbitron 
model; and, the orbitron theory is no longer considered as the 
correct description of the physics behind the observations re- 
ported by Alexeff and Dyer. 

Somewhat similar to the aforementioned pen sized coax- 
ial configuration investigated by Alexeff and Dyer is the mi- 
crohollow cathode discharges (MHCDs) configuration, which 
was first introduced by Schoenbach et al. 13 However, unlike 
the coaxial configuration, where the self-pulsing frequen- 
cies of ions typically lie in the gigahertz ranges, the self- 
pulsing of ions in the DC glow discharges for typical MHCDs 
configurations lie only in some tens of kilohertz frequency 
ranges . 21 i 51 i 52 Why? The answer to such discrepancy lies in 
the differences in the gas pressures applied in two configura- 
tions. In the high frequency design by Alexeff and Dyer, the 
pen sized coaxial device is filled with gas at a very low gas 
pressure of 0.1 mTorr whereas, in typical MHCDs configura- 
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tions, the device is maintained at gas pressures on the order of 
tens (or hundreds) of torrs. For instance, in the MHCDs ex- 
periment by Aubert et al. , 51 i 52 the gas pressure ranged from 40 
to 200 Torr. He et al£- worked with somewhat lower gas pres- 
sure of ~ 1 Torr (i.e., 133 Pa) for their MHCDs experiment, but 
compared to the gas pressure of . 1 mTorr used by Alexeff and 
Dyer in their pen sized coaxial configuration, this is still larger 
by a factor of ten thousand. Consequently, in the aforemen- 
tioned MHCDs experiments, the DC glow discharges involve 
much weaker static electric fields due to larger screening ef- 
fects. When a gas filled medium is applied with an external 
static electric field, the polarization process sets in and such 
screening effect weakens the net electric field strength in the 
region filled with gas. The degree of such screening process 
grows with gas pressure. Consequently, in typical MHCDs 
experiments, where much higher gas pressures are involved, 
relatively low self-pulsing frequencies are observed as a re- 
sult of larger screening effects for the applied static electric 
field. Contrary to this, in the aforementioned pen sized coax- 
ial configuration investigated by Alexeff and Dyer, the exter- 
nally applied static electric field is only weakly screened due 
to the fact that the device is maintained at a very low gas pres- 
sures. Consequently, the ions in such device can oscillate at 
very high frequencies. 

Besides this screening effect which acts to weaken the ex- 
ternally applied static electric field, the presence of higher gas 
pressure increases the effective mass of an oscillating super- 
particle. The superparticle concept has been previously intro- 
duced in the discussion of self-oscillations involving a plasma 
fireball. The dynamics of plasma involves collective motions 
of its constituent atoms. In the framework of this alternative 
theory 27,29 elaborated in this paper, the self-pulsing dynamics 
in both MHCDs experiments and the pen sized coaxial con- 
figuration investigated by Alexeff and Dyer involves the con- 
cept of oscillating superparticles. The effective masses of such 
superparticles have an explicit dependence on gas pressures. 
In general, the superparticle associated with the higher gas 
pressure environment has a larger mass than the one associ- 
ated with the lower gas pressure environment. Consequently, 
the superparticle in a typical MHCDs experiment has much 
larger mass than the one in the pen sized coaxial configuration 
investigated by Alexeff and Dyer. This also explains why self- 
pulsing frequencies are much lower in MHCDs experiments 
compared to the device investigated by Alexeff and Dyer. 



V. CONCLUDING REMARKS 

In summary, the mechanism behind self-sustained oscilla- 
tions of an ion in the DC glow discharges has been briefly dis- 
cussed in the framework of interaction between an ion and sur- 
face charges that it induces at the bounding electrodes. Such 
alternative description provides an elegant explanation to the 
formation of plasma fireballs in the laboratory plasma. It 
has been found that oscillation frequencies also increase with 
ion's surface charge density, but at the rate which is signifi- 
cantly slower than it does with electric field strength. Such 
result supports the conclusions by Fox 3 and Bosan et ai., 36 



where they reported of oscillation frequencies being not too 
dependent on the type of ions used in the discharge. Self- 
sustained oscillations in the DC glow discharges can be ex- 
pected to have high thermal stability due to ion's kinetic en- 
ergies that are much larger than the average room tempera- 
ture thermal energies. It is well known that self-sustained 
oscillations in the DC glow discharges occur only in the 
negative resistance region of the voltage-current characteris- 
tic curve. Experimentally, such region is characterized by 
Vthi < Vt < Vtfi2, wnere Vthi and V t h2 are some threshold volt- 
ages. Such observation is quite elegantly explained from the 
solutions of Eq. (0, where the physical solutions are only 
found for Vjhi < Vr < Vth2- Presented mechanism also cor- 
rectly describes the self-sustained oscillations of ions in dusty 
plasmas. To demonstrate this, Eq. (0 has been applied to 
correctly predict the frequency of dust particle oscillations 
in the dusty plasmas experiment by Arnas et a/.— Such re- 
sult demonstrates that self-sustained oscillations in dusty plas- 
mas and DC glow discharges involve common physical pro- 
cesses. This is the first theory to successfully explain the self- 
sustained oscillations phenomena in both dusty plasmas and 
the DC glow corona physics. 



APPENDIX 

A. Sheath potential 

It is worthwhile to compare the potential energy well of 
Fig. 0with the empirical potential well introduced by Tomme 
et al^ to describe the potential in the plasma sheath. The 
plasma sheath is an empty space residing between the plasma 
and the confining electrodes. One such plasma sheath near 
the anode is illustrated in Fig. a). For the reason that elec- 
trons and ions move at different velocities due to the differ- 
ence in their masses, plasmas are never completely neutral at 
any instant. Consequently, a plasma confined between DC 
voltage biased electrodes effectively behaves as if it were a 
super large charged-particle. For brevity, such "super large 
charged-particle" shall be simply referred to as a "plasma 
ball" throughout the discussion here. In fact, the single ion- 
ized nanoparticle case illustrated in Fig. 0b) is the limit in 
which the plasma in Fig. 0a) reduces down to contain just 
one ionized nanoparticle. That said, just as single ionized 
nanoparticle self-oscillates when confined by DC voltage bi- 
ased electrodes, the plasma ball also self-oscillates between 
the anode and cathode electrodes. 9 However, due to its large 
mass, the plasma ball oscillates at much lower frequencies 
compared to the nanoparticle counterparts. For instance, as- 
suming the gap h between the plate electrodes is in the order 
of sub-millimeters or so in Fig. 0a), the macroscopic plasma 
confined between such electrodes can contain very large num- 
ber, i.e., say, millions or more, of ionized nanoparticles or 
atoms depending on the gas pressure. In general, for the study 
of self-sustained oscillations, the plasma ball illustrated in Fig. 
0a) can be effectively modeled using an ionized superparticle, 
where the terminology superparticle refers to a particle which 
is extremely large and enormously heavier compared to the 
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ionized nanoparticle (or atom) illustrated in Fig. EJb). In such 
model, the potential in the plasma sheath is represented by 
the one presented in this article. It is quite remarkable that 
the sheath potential introduced empirically by Tomme et al. 
closely resembles the potential described in this work. 



(b) 



(a) 








K 3 


Sheath 


J z d 



Sheath 




Ionized nanoparticle 



E„ 



0.5h 



Figure 5: (Color online) (a) Illustration of the plasma and the sheath. 
The plasma ball effectively behaves as one very large charged super- 
particle. In the illustration, the plasma has been deliberately drawn as 
a sphere to emphasize a plasma ball. In reality, however, the plasma 
can be in any shape depending on the geometry of confining elec- 
trodes, (b) The single ionized particle configuration considered in 
this article. In (a) and (b), only the sheath near the anode is shown 
for illustration. 



by 02 = Ql (47ib 2 ) or 02 = — 5.35^lC/m 2 . Arnas et al. does 
not provide any information regarding the dielectric constant 
K2 for their glass microsphere. Therefore, K2 = 6.5 has been 
chosen for the dielectric constant of glass microsphere, which 
value is typical of glass microspheres. For convenience, it 
shall be assumed that the medium in which the glass micro- 
sphere oscillates is a vacuum, i.e., K3 = 1. The parameters in 
the experiment by Arnas et al. is summarized here for refer- 
ence: 



K 2 



■ 6.5, fc 3 = 1 , h 
rji =0C/m 2 , 
m rs 



5mm. E p = 5.15kV/m, 



o 2 = -5.35 juC/m 2 , 
15.1ng, a = 0m, b = 32jxm, 



(10) 



where the mass is in nanograms, i.e.. 
15.1 ng. The initial conditions, 



1.51 x l(r n kg = 



Za(0)= 4.4mm and 



dzj (0) = mm 



dt 



, (ID 



have been chosen purely out of convenience. With Eqs. ([Tot 
and (fTTT i. the equation of motion defined in Eq. © can be 
solved numerically via Runge-Kutta method to obtain the re- 
sults illustrated in Fig. |4] One can readily verify that the the- 
ory agrees with the experiment^ 



B. Parameters in experiment by Arnas et al. 

Arnas et al. 40 reported that a dust particle made of hollow 
glass microsphere of radius r 32 jim with mass density p m m 
110kg/m 3 and carrying a total charge of Q m —4.3 x I0 5 q e , 
where q e = 1.602 x 10~ 19 C, self-oscillates at frequency of 
V sc ~ 17 Hz inside the plasma sheath. The inner environ- 
ment of the plasma sheath, in many respects, is similar to the 
environment in an empty space between the DC voltage bi- 
ased plane-parallel conductors. In that regard, the undamped 
self-sustained oscillations of charged glass microsphere in the 
experiment by Arnas et al. should be theoretically describ- 
able using the model presented in this article, where a core- 
shell structured charged particle is confined between the DC 
voltage biased plane-parallel conductors. The glass micro- 
sphere used by Arnas et al. is not core-shell structured, of 
course. Such dielectric particle is obtained in the limit the 
radius of conductive core portion of the core-shell structured 
particle goes to zero, i.e., a = 0. Consequently, Ci also van- 
ishes in that limit. That said, the glass microsphere of radius 
b = 32jUm with mass density p„, » 1 10kg/m 3 has total mass 
of m = (4/3) 7ib 3 p m or m « 1.51 X 10~ n kg. Although Arnas 
et al. specifically uses the word "hollow glass microsphere" 
in their report, they do not provide any physical details of the 
particle other than its outer radius. Thus, in the aforemen- 
tioned calculation of the particle's mass, I have assumed a 
solid glass microsphere. The glass microsphere carries a to- 
tal charge of Q « -4.3 x I0 5 q e , where q e = 1.602 x 10~ 19 C. 
The surface charge density at the radius r = b is hence given 



C. Comparison to the potential energy of negatively charged 
glass microsphere measured by Arnas et al. 

According to the model elaborated here^ 9 - the positively 
charged particle confined between a DC voltage biased plane- 
parallel conductors results in the formation of potential energy 
well in vicinity of the anode whereas a negatively charged par- 
ticle results in the formation of potential energy well in vicin- 
ity of the cathode. Arnas et al. i5 have experimentally verified 
such potential energy well for the case of negatively charged 
particle in the plasma sheath near the cathode. As explained 
previously, the problem of charged particle inside the plasma 
sheath can be effectively modeled by a charged particle con- 
fined by the DC voltage biased plane-parallel conductors. The 
potential energy function U (z c i) of Eq. has been plotted for 
the following parameter values: 

!K"2 = 6.5, K3 = 1 , /i = 5mm, E p = 3.15kV/m, 
CFj =0C/m 2 , a 2 = -5.795 £iC/m 2 , (12) 
m «4.91 x 10~ 12 kg, a = 0m, b = 22\im. 

The result is shown in Fig. [6] where it shows the formation of 
potential energy well in vicinity of the cathode. One notices 
that the order of magnitude for the potential energy well is 
comparable to the experiment by Arnas et al. The minor dis- 
crepancies in the potential energy well of Fig. [6] and the one 
measured by Amas et alM can be attributed to the fact that the 
electric field is not constant in the plasma sheath whereas, be- 
tween the DC voltage biased plane-parallel conductors, elec- 
tric field is a constant. 
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Figure 6: The potential energy of Eq. ((4} is plotted for parameter 
values defined in Eq. d 1 2b . In the plot, the anode is located at Zd = 
0mm and the cathode is located at z d = 5 mm. This result should be 



compared with the potential energy measured by Arnas et a/— in 
their dusty plasmas experiment. One can verify that two results are 
remarkably similar in potential energy profile as well as in the order 
of magnitudes. 



D. Physical and unphysical oscillatory solutions 

The oscillatory solutions obtained for E p < 0. 1 GV /m and 
E p > 4.9 GV/m are unphysical in Fig. |3d). To demon- 
strate this, Eq. (0 is numerically solved and plotted for 
E p < 0.1 GV/m and E p > 4.9 GV/m using the following ini- 
tial conditions and parameter values: 



K-2=°°, K3 = l, /i = 100nm, V L = 0V, 
(7i «3.19C/m 2 , a 2 = 0C/m 2 , 
m w 8.95 x l(T 24 kg, b = a = 1 nm, 
Zd(0)=0.5nm, Zd(0)=0nm/s, i d = dzdldt. 



(13) 



To show that E p < 0.1 GV/m yields unphysical solutions, 
Zd (t) is plotted for E p = 0.1GV /m and E p = 0.09 GV /m. The 
results are shown in Fig. [7] where the surface of anode is at 
Zd = Onm, the surface of cathode is at Zd = lOOnm, and the 
midway between the plates is at Zd = 50 nm. It can be clearly 
seen that for E p = 0.09 GV / m, positively ionized particle peri- 
odically crosses the midway between the two electrodes. Such 
case is unphysical because, for positively ionized particles, 
oscillations can only exist for Zd < h/2 — Now, it can be 
verified that any E p less than E p = 0.09 GV/m yields such 
unphysical solutions. 

To show that E p > 4.9 GV/m yields unphysical solutions, 
z d (t) is plotted for E p = 4.9GV/m and E p = 5.0GV/m. The 
results are shown in Fig. HJa), where the plot has been en- 
larged for a view near Zd = 0pm in Fig. [Sb)- It can be 
clearly seen that for E p = 5.0GV/m, positively ionized par- 
ticle periodically penetrates into the surface of anode, which 
is unphysical. Now, it can be verified that any E p larger than 
E p = 5.0 GV/m yields such unphysical solutions. 

For the initial conditions and parameter values defined in 
Eq. (jT3j, physical solutions for the oscillatory motion of pos- 
itively ionized particle are only found for E p satisfying the 
condition given by 0.1 GV/m < E p < 4.9GV/m, where the 



Figure 7: (Color online) Plot of Zd(t) for E p = O.lGV/m and 
E p = 0.09GV/m. The surface of anode is at Zd = Onm, the surface 
of cathode is at Zd = lOOnm, and the midway between the two elec- 
trodes is at Zd = 50nm. 
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Figure 8: (Color online) (a) Plot of Zd (0 for E p = 4.9 GV/m and 
Ep = 5.0 GV/m. (b) The plot in (a) has been enlarged for a view 
near z d = 0pm. The surface of anode is at Zd = 0pm. 



upper and lower bounds of E p are rough estimates based on 
the grounds of the discussed physicality. This result explains 
why oscillations suddenly appear at certain "initial" threshold 
DC bias voltage and disappear suddenly when bias voltage 
goes beyond certain larger "final" threshold voltage in exper- 
iments. 



10 



E. Weakly ionized particle confined between anode and 
cathode plates separated by gap of h = 1 mm 

The oscillation frequency plot of Fig. |5Jd) has been ob- 
tained for an highly ionized particle which is confined be- 
tween two plate electrodes with very small separation gap of 
h = lOOnm. Here, the same calculation is done for a weakly 
ionized particle confined between plate electrodes with mi- 
croscopically very large separation gap of h = 1 mm. Ac- 
cording to Chen et al.,— a spherical argon cluster of radius 
r = a = 4.7 nm contains roughly ~ 9600 argon atoms at back- 
ing pressure of 50 bars. Since the particle is assumed to be 
weakly ionized, I shall assume that N — 20 electrons are re- 
moved from it. This corresponds to surface charge density 
of ai = Nq e / (4-Tta 2 ) or ai « 1.1543 x 10- 2 C/m 2 , where 
q e = 1 .602 x 10~ 19 C is the fundamental charge unit. Neglect- 
ing the masses of missing N electrons, the particle has a total 
mass of m m 6.365 x 10~ 22 kg. For the calculation of particle's 
mass, it has been assumed that single argon atom has mass 
of 6.63 x 10~ 26 kg. That said, the following initial conditions 
and parameter values are used for the calculation of particle's 
oscillation frequency as function of electric field strength: 

K2 = °°, K3 = 1 , h = 1mm, Vl = V, 

cti w 11.543mC/m 2 , a 2 =0C/m 2 , 
m « 6.365 x 10~ 22 kg, b = a = 4.7nm, { ' 

z^(0) = 10jUm, td{0) =0Aim/s, z d =dz d /dt. 

Using these, the equation of motion defined in Eq. © is 
solved numerically via Runge-Kutta method for the oscilla- 
tion frequency as function of electric field strength; and, the 
result is shown in Fig. [9] 

This result should be compared with the one provided by 
Akishev et al., 25 where they have plotted both experimen- 
tal and calculated period r osc of self-sustained oscillations 
against the average corona current at different pressures of 
nitrogen. Fox 3 and Bosan et al.— have shown that oscilla- 
tion frequencies are not very dependent on the type of ionized 
particles used in the discharge. That said, the result obtained 
by Akishev et al. for nitrogen gas can be compared with the 
calculation done here for ionized spherical argon cluster. Re- 
calling that v osc = 1 /Tose and the average corona current in 
the electrode increases with electric field strength, it can be 
readily convinced that the profile of result shown in Fig. |9]is 
consistent with the result obtained by Akishev et al. 

Although the profile of oscillation frequency dependence 
on electric field strength (or corona current in the electrode) 
is comparable in both results, the measurement by Akishev et 
al. shows much lower oscillation frequencies for given range 
of electric field strengths. Why? Such discrepancy arises 
from the fact that in the experiment by Akishev et al., the 
self-oscillating object is a charged plasma ball, i.e., charged 
superparticle, whereas in the calculation of Fig. [9] the self- 
oscillating object is a charged nanoparticle. As already dis- 
cussed in section A, a macroscopic plasma ball contains very 
large number of ionized nanoparticles (or atoms) depending 
on the gas pressure. This makes macroscopic plasma ball ef- 



fectively a single charged superparticle with very large mass. 
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Figure 9: (a) The frequency of oscillations is plotted against electric 
field strength for a weakly ionized particle with initial conditions and 
parameters values defined in Eq. l!14t , (b) The plot in (a) is enlarged 
near zero. Unlike the case illustrated in Fig. Hid), where ionized 
particle is confined between electrodes with very small gap, the re- 
sult here involves very large number of discretized time steps for 
the Runge-Kutta routine. For that reason, the oscillation frequency 
has been plotted from E p = 0.01 kV/m to E p = 3.5kV/m in (a). The 
threshold electric field strengths, i.e., Ep tbl an d £p,th2 > 316 not shown 
in the figure. 



It has been illustrated in Fig. [4] that an ionized particle with 
smaller mass requires weaker electric field strength compared 
to the one with larger mass to oscillate at the same frequency. 
Mathematically, such characteristic arises from the fact that 
Eq. © has mass dependence in the denominator. Conse- 
quently, the plasma ball in the experiment by Akishev et al., 
which effectively behaves as a single charged superparticle, 
oscillates at much lower frequencies for the given range of 
electric field strengths compared to the configuration used in 
Fig. [9] where the mass of an ionized particle is much smaller 
than the effective mass of a plasma ball. In principle, once the 
mass and the effective total charge information of the plasma 
ball, i.e., superparticle, are provided, the oscillation frequency 
dependence on corona current in the electrode (or electric field 
strength) measured by Akishev et al. can be reproduced by the 
presented theory. 
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